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Previous studies have shown that noise can induce coherence resonance in some nonlinear dynami- 
cal systems close to a bifurcation of a periodic motion, such as in excitable systems. We demonstrate 
that coherence resonance can be observed in systems close to a blowout bifurcation. It is shown 
that for dynamical systems with an invariant subspace in which there is a phase-coherent chaotic 
attractor, the interplay among the oscillation of local transverse stability, noise and nonlinearity can 
lead to coherence resonance phenomenon. The mechanism of coherence resonance in this type of 
system is different from that in previously studied systems. 

PACS number(s): 05.40.-a, 05.45.-a 
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I. INTRODUCTION 

The behavior of nonlinear dynamical systems sub- 
jected to noise has been an interesting subject of recent 
investigation. A great deal of work has been devoted 
to stochastic resonance jjj, where an optimal amount of 
additive noise can generate the maximal response of the 
system to a weak external periodic or aperiodic signal. 

Coherent motion can be induced purely by noise in 
some dynamical systems without an external signal. In 
nonlinear dynamical systems near the onset of bifurca- 
tions of periodic orbits, the periodicity is visible even 
before the bifurcation actually occurs if there is noise 
present, a phenomenon called noisy precursor of the bi- 
furcation H . Similar phenomenon has also been observed 
in excitable systems ||, such as in various neural mod- 
els and laser system fjj, in the fixed point regime 
close to a saddle-node bifurcation of a periodic orbit. A 
feature common to the systems close to the bifurcation 
of a periodic orbit is that during the relaxation to sta- 
ble orbit below the bifurcation, the transient possesses 
the periodicity above the bifurcation, and the effect of 
external noise is to continually kick the system off of the 
stable orbit, so that the transient behavior displays co- 
herent motion. 

More interestingly, recent investigations have shown 
that the coherence of the noise-induced motion achieves 
a maximum at an optimal noise intensity For ex- 

ample in an excitable system, when the system is kicked 
away from the fixed state to overcome a certain thresh- 
old, it will come back to the fixed point only after a large 
excursion (noise-induced limit cycle, or spike in neural 
systems). When noise is weak, the system is rarely ex- 
cited, and the motion is quite irregular. Increasing noise 
kicks the system over the threshold more often, and the 
system fires more and more spikes. The interspike inter- 
val becomes the most regular at an optimal noise level. 
After that, too high level of noise distorts greatly the near 
limit cycle, rendering the motion irregular again. Simi- 



lar to the conventional stochastic resonance, this phe- 
nomenon of resonance without an external signal is called 
coherence resonance (CR). In a recent paper, Ohira and 
Sato showed that a simple two state model with time 
delay can display CR |L3]]. Indeed, the time delay in- 
troduces an intrinsic periodic oscillation into the system: 
if noise induces a spike at a certain moment, another 
spike is mostly expected after the time delay. More re- 
cent work has also shown noise-enhanced synchroniza- 
tion ||-|ll| and array-enhanced CR G3| in coupled or 
extended excitable systems. 

It is interesting and practically meaningful to see 
whether CR exists in other type of system. In this pa- 
per, we demonstrate CR in noisy chaotic dynamical sys- 
tems close to a blowout bifurcation |l4|,|l5| which occurs 
in dynamical systems with an invariant subspace S in 
the phase space. The transverse stability of the sub- 
space is determined by the motion within the subspace. 
When the largest transverse Lyapunov exponent A is neg- 
ative, S is stable; while unstable when A is positive, 
and the critical point of the transverse stability is the 
blowout bifurcation point. However, the local stability 
of the subspace may fluctuate greatly when the motion 
within the subspace is chaotic. The finite time Lyapunov 
exponent A^ measuring expansion or contraction of a 
transverse perturbation during a period of finite time T 
may oscillate greatly around the average value A. This 
fluctuation of local stability can lead to interesting and 
unusual behaviors, such as bubbling jlj], on-off inter- 
mittency Jl5[ and additional complexity in the system 
by unstable dimension variability jfujj . When there is 
a quasiperiodic torus in S, the system may undergo a 
transition to strange nonchaotic attractors via blowout 
bifurcation 11171 - 



1 



II. SYSTEM AND RESULTS 

CR observed in this type of system is the most appre- 
ciable when the motions within the subspace are oscilla- 
tions with property of being phase- wise, although it may 
be chaotic in the amplitude. A typical example of this 
type of motion is the Rossler chaotic attractor [Fig. 1] . 
In this work, we focus on the case where there is such a 
motion in the subspace, and study the system behavior 
in the presence of noise. Since the noise prevents the dy- 
namics from approaching the subspace indefinitely the 
system may be repelled far away from the subspace due 
to the local instability. Previous investigations on the ef- 
fects of noise in this type of system focused on the change 
of the universal behavior of the laminar phase distribu- 
tion |Q of the blowout motion in the context of on-off 
intermittency. Here we focus on the coherence of the 
blowout motion. Our results will show that, with the 
increase of the noise level, the output (e.g. the distance 
from the subspace) displays increasing coherence till too 
high level of noise dominates the dynamics and destroys 
the coherence, exhibiting typical CR phenomenon. To 
illustrate our findings, we consider the case where there 
is a Rossler chaotic attractor in the invariant subspace, 



x\ = a(-x 2 - x 3 ), 
x 2 — a{xi + ax 2 ), 
x 3 = a(0A+ (x 1 -8.5)0:3), 
y = [b(xi - xi) + c] sin y - y - 



(1) 

(2) 
(3) 
(4) 



where a controls the time scale of the Rossler system, 
and x\ is the time average of x\. Clearly y — defines 
the invariant subspace in the absence of noise £ which is 
a Gaussian white noise with (£,(t)t;(t')) = 5(t — t'). The 
introduction of x\ in Eq. (4) is for the convenience of 
discussion, because the largest transverse Lyapunov ex- 
ponent 



A = lim - 

t—>oo t 



[b(xi — xi) + c — l]dr 



1. 



(5) 



and the transverse stability is only controlled by the pa- 
rameter c. Eqs. (1-4) are modified version of a physical 
model of superconducting quantum interference device 
(SQUID) Jl9| . Similar model with a torus in the subspace 
has been studied in the context of nonchaotic strange at- 
tractors In fact, the specific form of the nonlinearity 
in Eq. (4) is of no importance for the CR phenomenon. 
It serves to keep the system bounded. 

We note that for parameter a = 0.15, the Rossler 
system possesses a chaotic attractor with some degree 
of phase coherence [see Fig. 1(a)]: the cycling time Tr 
has a rather sharp distribution, or equivalently, there 
is a pronounced peak at frequency luq on the broad- 
band spectrum of the chaotic signal. Phase synchroniza- 
tion [^,|T) and lag synchronization |^2j occur in coupled 
Rossler systems due to this phase coherence. For the 
parameter c near the blowout bifurcation point c = 1.0, 



the local stability undergoes large chaotic oscillations. 
The subspace is only temporally attracting in about half 
period of the loops when b{x\ — 5i) + A < 0, while 
is temporally repelling in the other half period when 
b{x\ — x\) + A > 0. In the following, we study sys- 
tem behavior in the presence of noise both below and 
above the blowout bifurcation c = 1.0, with parameters 
a = 0.1,a = 0.15,6=\/0l)2. 
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FIG. 1. Chaotic attractors of Rossler system with different 
parameters (a) a=0.15 and (b) a=0.25. 



A. Below the blowout bifurcation: c < 1 

Typical behaviors of the output |y| in the presence of 
noise with different levels are shown in Fig. 2 for c = 0.9. 
When the noise is very weak, the system stays very 
closely to the invariant space y = 0, giving no appre- 
ciable output [Fig. 2(a)]. If the noise is larger than a 
certain level, the system begins to produce large output, 
but the signal is quite irregular [Fig. 2(b)]. At a certain 
range of noise levels, the output become very regular; it 
is almost periodic, with the amplitude fluctuating only 
slightly [Fig. 2(c)]. However, when noise goes to even 
higher levels, it begins to deform the near periodic sig- 
nal, and degrades the regularity [Fig. 2(d)]. The response 
property of the system to additive noise exhibits typical 
feature of CR: the system motion becomes the most co- 
herent at an optimal noise level. 

To characterize the degree of coherence, we compute 
the spectra of the output signals \y\. It has a pronounced 
peak at the same frequency c^o as the driving signal x\. 
Fig. 3(a) shows the peak height p m at ujq as a function 
of the noise level D\ — log 10 u. Below a certain value of 
Di, p m increases exponentially and quickly; then it comes 
to a slow increasing region covering many orders of the 
noise level. p m reaches a maximal value and decreases at 
higher noise levels. 

In the following, let us look into the mechanism of CR 
in the system. In the absence of noise, the subspace is 
stable, and from a random initial condition, \y\ decreases 
exponentially on average as \y\ ~ exp(At), but with a lo- 
cal chaotic oscillation. The additive noise, however, pre- 
vents the system from approaching the subspace much 
deeper than the noise level, and the trajectory can be 
repelled away from the noise level to produce relatively 
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large output. Thus the noise sets a reflection boundary 
to the transverse motion, as seen in the corresponding 
semilog plots of the output in the right panel of Fig. 2. 
Due to this boundary, the transverse motion is attracted 
to and repelled away from the noise level alternately when 
the motion within the subspace comes into the local sta- 
ble region b(xi — x\) + A < and local instable region 
b(xi — xi ) + A > 0. The intrinsic coherence of the chaotic 
motion within the subspace is thus manifested by the 
noise-induced blowout motion. This boundary effect of 
noise is quite different from that in excitable systems and 
the systems with delay, where noise acts to kick the sys- 
tem over a threshold. 



\y\ w exp[A(cosuj t + 1) + D] 



(7) 
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FIG. 2. Typical behavior of noise induced blowout motion 
(left panel) at different noise levels a. The right panel is the 
corresponding semilog plots of the same quantity. For the 
purpose of clear illustration of the boundary moving with the 
noise level, we use the same scale in the right panel. The 
dotted lines in the right panel indicate the noise level, (a) 
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(b) a = 10", (c) a = 10 , and (d) a = 10" 



The calculation of p m is not easy from solving the 
complete stochastic system in Eq. (4) which is driven 
multiplicatively by chaotic signal x\. We choose to esti- 
mate it with more intuitive but quite accurate approxi- 
mations based on the above observation that the effect 
of noise can be modeled by a reflection boundary. For 
weak enough noise, \y\ keeps small, and above the noise 
level \y\ 3> a, the dynamics is governed by the linear 
approximation y — [b(xi — x±) + A]y, which gives 



exp 



[b(xi -xi) + A]dr 



(6) 



Now it is not surprising that the spectra of \y\ have a 
peak at the same frequency ujq as x\. Note the intrin- 
sic coherence of the chaotic signal x±, we can roughly 
approximate x\ by a periodic signal with frequency ojq. 
Take into account the boundary of the noise, the right 
panel of Fig. 2 suggests that we may approximate 
f Q [b(xi — xi) + A]dr ss A(cosuJot + 1) + D, and it fol- 
lows from Eq. [j] that 



with a proper shift of the time origin. The amplitude 
A = [b(x m — x\) + Air]/u>o, where x m is the average 
of the maxima of x\. 2 A is the average of the maxi- 
mal magnitude that \y\ can depart from the boundary 
at the noise level D = In a during a cycle of x\ . With 
the above system parameters, we numerically estimate 
x m = 11.42, x\ = 0.1324 and wo = 1.035a, thus giv- 
ing A = 12.40. As the noise level increases, the bound- 
ary is moving to a higher order [Fig. 2]. If —D < 2 A, 
the maximal value of \y\ comes to the order of unit, and 
the system begins to produce relatively large output, and 
the nonlinearity begins to set in. However, due to chaotic 
fluctuation of the amplitude of x\, the maxima of \y\ may 
not be as large during those small cycles of xi[Fig. 2(b)], 
and the behavior is quite irregular, as is familiar in the 
context of on-off intermittency. With the boundary going 
to even higher order, the system can also produce quite 
large output \y\ during those small cycles of x\ and the 
amplitude is confined to some saturated values by the 
nonlinearity [Fig. 2(c)]. The system now performs very 
coherently. In this nonlinear regime, Eq. p based on lin- 
ear approximation is no longer valid. However, the fact 
that \y\ always begins to rise from the noise level when x\ 
comes into local instable region b{x\ — x\) + A > shows 
that the blowout motion still possesses the frequency lvq 
in the strong nonlinear regime. In addition, Fig. 2(c) 
suggests that we can still approximate In \ y\ by a peri- 
odic function whose amplitude is between the boundary 
of noise level and the confinement of the nonlinearity, i.e., 



\y\ 



exp[£?(cosojo£ + 1) + D\. 



(8) 



Now the amplitude B — (D m — D)/2, where D m = 
lnmax(|?/|) « 0.8 with the above system parameters. 

With the approximations in Eq. ^ and Eq. ^, we can 
estimate the peak height at wo as a function of D in dif- 
ferent dynamical regimes as 



2tt 



■tt/ujo 

I y | cos coot dt 

F 2 (A)e 2A e 2D , -D > 2A, 
F 2 (B)e 2B e 2D , -D < 2A. 



(9) 



Here, F(A) = J^cost exp(Acost) dt. The analyt- 
ical result explains the exponential increase of the peak 
height in the linear regime observed in numerical simu- 
lation [Fig. 3(a), cycles]. With A = 12.80, Eq. | fits the 
exponential region very well [Fig. 3(a), solid line]. Note 
that the fitting parameter A = 12.80 is quite close to 
A = 12 AO estimated from the system parameters, show- 
ing that the above simple approximations give a good 
account for the system behavior. The analysis also re- 
produces qualitatively the result in the nonlinear regime. 
The discrepancy in the crossover region (— D ~ 2A) is 
due to the fact that in this region the maximal value of 
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\y\ is not saturated and D m = 0.8 used in fitting over- 
estimates the maximal values of \y\ in this region. One 
should note that in the above analysis, the specific form 
of nonlinearity is of no importance, and only the confine- 
ment property of the nonlinearity is employed, indicating 
that the phenomenon is universal in this type of systems. 
Now we understand that CR occurs in the system due to 
the interplay among the chaotic while somewhat coherent 
oscillation of the local stability of the subspace, the con- 
finement of the nonlinearity to the transverse motion and 
the boundary effects of the additive noise. This mecha- 
nism is quite different from that in previously investi- 
gated systems [p|~|l3|]. 
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FIG. 3. Illustration of coherence resonance below the 
blowout bifurcation, (a) Peak height p m as a function of noise 
level D\ = log 10 a. The solid line is the analytical estimation 
with A = 12.80. (b) Coherence factor /3 as a function of noise 
level. The solid line is the result of the chaotic signal xi. 

Based on the above understanding of the behavior, we 
can characterize CR by another quantity, the "coherence 
factor" defined by the relative fluctuation of the output 
amplitude, e.g. the time average of the amplitude A y 
divided by its standard deviation 

P=(A y )/yJvzr(A y ). (10) 

For relatively weak noise, \y\ 3> a most of the time, and 
\y\ is smooth and A y is well defined. For relatively large 
noise, noise-induced short time fluctuation of \y\ can be 
rather strong. In numerical simulations A y is defined as 
follows: firstly a smooth series \y\ m is obtained from the 
noisy function \y\ using a moving average method, then 
A y is taken as the value of \y\ at the moment when \y\ m 
is maximal. This definition captures the noise-induced 
short time fluctuations of \y\ at large noise levels. The re- 
sult is shown in Fig. 3(b) (cycles). (3 reaches a maximum 
and decreases quickly when too high level of noise begins 
to dominate the fluctuation of the amplitude. In linear 
and weak nonlinear region (Z?i < —7.5), the chaotic fluc- 
tuation of the amplitude of x\ is augmented by the expo- 
nential relationship between \y\ and x\, as seen in Eq (|6|), 
so that the coherence degree of \y\ is lower than that of 
xi, as can be seen by the comparison of the coherence 
factors of \y\ (cycles) and x\ (solid line). It is very in- 
teresting to see that in the nonlinearity dominant region, 
the coherence of the noise-induced motion is much higher 
than the intrinsic coherence of x\ in a wide range of the 



noise level, because the confinement of the nonlinearity 
smoothes the fluctuation of the amplitude of \y\. Thus 
the combination of the noise and nonlinearity enhances 
greatly the intrinsic coherence. Such pronounced CR 
phenomenon is able to be observed when other chaotic at- 
tractors possessing similar intrinsic coherence are within 
the subspace, such as the electronic circuit in Ref. [23]] 
and the hybrid laser system in Ref. [^4j or the ecological 
system in Ref. |^5| . It is also clear that CR will occur for 
periodic and quasiperiodic motion within the subspace, 
i.e., the peaks in the spectra possess a maximal value at 
an optimal noise level. 

In the Rossler system, the topology of the chaotic at- 
tractor changes if a is large than 0.21: there are large 
and small loops [see Fig. 1(b), a=0.25]. This topology is 
quite typical in many low dimensional chaotic systems. 
Now, due to the large fluctuation of the amplitude and 
returning time, there are no pronounced peaks in the 
broadband spectra of the chaotic signals. However, CR 
is still observable. Fig. 4 shows the coherence factor (3 for 
| y | and x\. Again we see the enhancement of the intrin- 
sic coherence by noise, albeit with less intensity. Similar 
behavior should be observable in general systems of this 
type. 
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FIG. 4. Coherence factor j3 as a function of the noise level 
for chaotic motion with funnel Rossler attractor in Fig. 1(b). 



B. Above the blowout bifurcation: c > 1 

Typical behavior of the system above the blowout bi- 
furcation is similar to that below the bifurcation. The 
difference is that, for c > 1, the subspace is transversely 
unstable, and blowout motion has already existed even 
without noise, and the confinement of the nonlinear- 
ity has already taken place. For c < 1, the dynam- 
ics can always access the noise level, however weak it 
may be; while for c > 1, the dynamics can only come 
to the subspace no closer than \y\ ~ exp(— 2 A), where 
A = [b(x m — x\) — Ait]/ll)q for the chaotic attractor in 
Fig. 1(a). As a result, weak noise with —D> 2 A will 
have no discernible effects on the system behavior. How- 
ever, stronger noise will act as a reflection boundary simi- 
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lar to the case below the blowout bifurcation, preventing 
the system from approaching the subspace to the clos- 
est level \y\ ~ exp(— 2A). The system behavior is now 
very similar to that in the nonlinear regime below the 
bifurcation. This analysis is demonstrated by numeri- 
cal simulations with c = 1.1 in Fig. 5. It is seen that 
both p m and f3 keep unchanged for very weak noise, and 
the coherence increases once the noise becomes effective. 
While p m increases monotonically till the noise domi- 
nates the dynamics and destroys the coherence of the 
blowout motion, (3 also exhibits another peak at rather 
weak noise level (a <~ 1CP 10 ). This behavior is related 
to the topology of the chaotic attractor in the subspace. 
As seen in Fig. 1(a), the attractor always flips to the 
smallest loops from the largest ones. For weak noise 
without significant effects, the maxima of \y\ associated 
with those small loops have relatively small values, which 
contributes mainly to the fluctuation of the amplitude of 
\y\. The broad distribution of the amplitude is clearly 
illustrated by the return map of the maxima of \y\ for 
g = 1CT 14 in Fig. 6. For a <~ 1CT 10 , the system has 
access to the noise level when x\ cycles along the largest 
loops. The reflecting property of the noise has the effect 
to increase greatly the maximal values of \y\ for those 
smallest loops following those largest ones, while only 
slightly for others. Most of the maxima are confined to 
a small neighborhood around A y = 2.0, as seen by the 
crosses in Fig. 6. This reduces the fluctuation of the 
amplitude of \y \ and enhances the coherence greatly. For 
intermediate noise level a ~ 10~ 7 , the fluctuation be- 
comes a little larger again when the maxima of |y| are 
pushed to slightly larger values by the reflection of the 
noise. Further increase of noise pushes more maxima of 
|y| to fluctuate slightly around a saturated value till too 
high level of noise destroys the coherence, resulting in an- 
other peak. The return map close to the peak is shown 
in Fig. 6 for a = 1CT 2 . 
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FIG. 5. Illustration of coherence resonance above the 
blowout bifurcation, (a) Peak height p m as a function of noise 
level D\ = log 10 a. Unlike Fig. 3(a), linear scale is used for 
p m here, (b) Coherence factor f3 as a function of noise level. 
The solid line is the result of the chaotic signal xi. 
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FIG. 6. Return map constructed by successive maxima of 



\y\ for three different levels of noise a = 10 
and a = 1CT 2 . 
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CR phenomenon in the system is most appreciable 
around the blowout bifurcation, i.e., in the on-off inter- 
mittency regime. In general, the maximal coherence is 
higher for the system above the bifurcation point. If c is 
far below the critical point, only large enough noise can 
induce blowout motion from the subspace, and the co- 
herence of the motion may have already destroyed by the 
noise. The peaks of p m and /3 become lower and narrower 
as c decreases from c = 1.0 and disappear when the sub- 
space becomes stable almost everywhere. For c far above 
the critical point, the system has only seldom close access 
to the subspace, and weak noise has no significant effects 
on the motion, while strong noise degrades the coherence 
of the motion away from the subspace. Typically, one 
observes that both p m and (3 keep unchanged for weak 
noise and begin to decrease for high enough noise when 
the subspace becomes unstable almost everywhere. 



III. CONCLUSION 

In summary, we have shown that the phenomenon of 
coherence resonance can be naturally observed in some 
nonlinear dynamical systems possessing an invariant sub- 
space, close to the blowout bifurcation where previous 
studies were often in very different context of on-off in- 
termittency. A link between these two distinct dynamical 
phenomena, CR and on-off intermittency, is the fluctua- 
tion of the local transverse stability of the subspace due 
to the oscillatory motion within the subspace. The noise- 
induced blowout motion manifests the coherent oscilla- 
tion within the subspace, while the confinement of the 
nonlincarity reduces the chaotic fluctuation of the am- 
plitude. The additive noise and the confinement of the 
nonlincarity to the transverse motion combine together 
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to manifest and enhance the intrinsic coherence of the 
motion within the subspace to the maximal degree. CR 
in this class of system with mechanism different from that 
in previously studied systems extends our understanding 
on nontrivial and positive effects of noise on nonlinear dy- 
namical systems, and could be physically and practically 
meaningful. 
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